^ 



Collective Flow signals the Quark Gluon Plasma * 

H. Stocker 

Institut fiir Theoretische Physik, Johann Wolfgang Goethe - Universitdt, 

Robert Mayer Str. 8-10, 60054 Frankfurt am Main, Germany and 

Frankfurt Institute for Advanced Studies (FIAS), 

r—s , Robert Mayer Str. 8-10, 60054 Frankfurt am Main, Germany 

o ■ 

»vj ' A critical discussion of the present status of the CERN experiments on charm dynamics and 

I hadron collective flow is given. We emphasize the importance of the flow excitation function from 

H , 1 to 50 A-GeV: here the hydrodynamic model has predicted the collapse of the vi-flow and of the 

, ^ ' W2-flow at ~ 10 A-GeV; at 40 A-GeV it has been recently observed by the NA49 collaboration. Since 

hadronic rescattering models predict much larger flow than observed at this energy we interpret this 

C^^ , observation as potential evidence for a first order phase transition at high baryon density pB- A 

detailed discussion of the collective flow as a barometer for the equation of state (EoS) of hot dense 
matter at RHIC follows. Here, hadronic rescattering models can explain < 30% of the observed 
^ , elliptic flow, V2, for pt > 2 GeV/c. This is interpreted as evidence for the production of superdense 

OO ■ matter at RHIC with initial pressure far above hadronic pressure, p > 1 GeV/fm'^. We suggest that 

the fluctuations in the flow, vi and V2, should be measured in future since ideal hydrodynamics 
predicts that they are larger than 50 % due to initial state fluctuations. Furthermore, the QGP 
coefficient of viscosity may be determined experimentally from the fluctuations observed. The 
connection of V2 to jet suppression is examined. It is proven experimentally that the collective flow 
is not faked by minijet fragmentation. Additionally, detailed transport studies show that the away- 

<;;l . side jet suppression can only partially (< 50%) be due to hadronic rescattering. We, finally, propose 

,.^ ' upgrades and second generation experiments at RHIC which inspect the first order phase transition 

in the fragmentation region, i.e. at ^b ~ 400 MeV (j/ « 4 — 5), where the collapse of the proton flow 
should be seen in analogy to the 40 A-GeV data. The study of Jet- Wake-riding potentials and Bow 
shocks - caused by jets in the QGP formed at RHIC - can give further information on the equation 
of state (EoS) and transport coefficients of the Quark Gluon Plasma (QGP). 
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I. OLD AND NEW OBSERVABLES FOR THE QGP PHASE TRANSITION 
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5t 1 Lattice QCD results .1, 2] (cf. Fig. ^| show a crossing, but no first order phase transition to the QGP for vanishing or 

small chemical potentials jjlb, i-e. at the conditions accessible at central rapidities at RHIC fuU energies. A first order 
phase transition does occur according to the QCD lattice calculations [l|, Q only at high baryochemical potentials or 
densities, i.e. at SIS-300 and lower SPS energies and in the fragmentation region of RHIC, y ~ 4 — 5 [^|^. The critical 
baryochemical potential is predicted 0,0| to be /i^ « 400 ±50 MeV and the critical temperature Tc ~ 150— 160 MeV. 
We do expect a phase transition also at finite strangeness. Predictions for the phase diagram of strongly interacting 
matter for realistic non- vanishing net strangeness are urgently needed to obtain a comprehensive picture of the QCD 
phase structure. Multi-Strangeness degrees of freedom are very promising probes for the properties of the dense and 
hot matter |^. The strangeness distillation process |a,l3 predicts dynamical de-admixture of s and s quarks, which 
yields unique signatures for QGP creation: high multistrange hyperon-/-matter production, strangelet formation and 
unusual antibaryon to baryon ratios ect. 

A. Thermodynamics in the T — /is plane 

A comparison of the thermodynamic parameters T and hb extracted from the UrQMD-transport model in the 
central overlap regime of Au-|-Au collisions ^ with the QCD predictions is shown in Fig 1, where the full dots with 
errorbars denote the 'experimental' chemical freeze-out parameters - determined from fits to the experimental yields - 
taken from Ref. 10]. The triangular and quadratic symbols (time-ordered in vertical sequence) stand for temperatures 
T and chemical potentials jjlb extracted from UrQMD transport calculations in central Au-f Au (Pb+Pb) collisions 
at RHIC (21.3 A-TeV), 160, 40 and 11 A-GeV _8j as a function of the reaction time (separated by 1 fm/c steps from 
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FIG. 1: The new phase diagram with the critical end point at fiB ~ 400 MeV, T « 160 MeV as predicted by Lattice QCD. In 
addition, the time evolution in the T — /^-plane of a central cell in UrQMD calculations 8j is depicted for different bombarding 

energies. Note, that the calculations indicate that bombarding energies Elab ~ 40 A-GeV are needed to probe a first order 
phase transition. At RHIC (see insert at the ^ib scale) this point is accessible in the fragmentation region only (taken from 9]). 



top to bottom). The open symbols denote nonequilibrium configurations and correspond to T parameters extracted 
from the transverse momentum distributions, whereas the full symbols denote configurations in approximate pressure 
equilibrium in longitudinal and transverse direction. 

During the nonequilibrium phase (open symbols) the transport calculations show much higher temperatures (or 
energy densities) than the 'experimental' chemical freeze-out configurations at all bombarding energies (> 11 A-GcV). 
These numbers are also higher than the critical point (circle) of (2+1) flavor - Lattice QCD calculations by the 
Bielefeld-Swansea-collaboration 2] (large open circle) and by the Wuppertal-Budapest-coUaboration 1] (the star 
shows earlier results from 1]). The energy density at Hc,Tc is in the order of w 1 GeV/fm'^ (or slightly below). 
At RHIC energies a cross-over is expected at midrapidity, when stepping down in temperature during the expansion 
phase of the 'hot fireball'. The baryon chemical potential jjlb for different rapidity intervals at RHIC energies has been 
obtained from a statistical model analysis by the BRAHMS Collaboration based on measured antihadron to hadron 
yield ratios [Ij. For midrapidity one finds ^b — 0, whereas for forward rapidities /is increases up to /is ~ 130 McV 
at y = 3. Thus, only extended forward rapidity measurement (y « 4 — 5) will allow to probe large /i^ at RHIC. The 
detectors at RHIC at present offer only a limited chemical potential range. To reach the first order phase transition 
region at midrapidity, the International Facility at GSI seems to be the right place to go. This situation changes 
at lower SPS (and top AGS) as well as at the future GSI SIS-300 energies: sufficiently large chemical potentials 
jjiB should allow for a first order phase transition 12] (to the right of the critical point in the (T^^b) plane). The 
transport calculations show high temperatures (high energy densities) in the very early phase of the collisions, only. 
Here, hadronic interactions are weak due to formation time effects and yield little pressure. Diquark, quark and gluon 
interactions should cure this problem. 



Open diarm and charmonia at SPS and RHIC 



Open charm and charmonium production at SPS and RHIC energies has been calculated within the HSD and 
UrQMD transport approaches ,13j using parametrizations for the elementary production channels including the 
charmed hadrons D, 5, D*,D*,Ds, Dg, Dt,Dt, J/*, *(25), X2c from NN and irN collisions. The latter parametriza- 
tions are fitted to PYTHIA calculations '14l| above y^ — 10 GeV and extrapolated to the individual thresholds, while 
the absolute strength of the cross sections is fixed bv the exoerimental data as described in Ref. ilSil: for orevious 
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FIG. 2: The J/*? suppression as a function of the transverse energy Et in Pb + Pb coUisions at 160 A-GeV. The sohd hue 
shows the HSD resuh within the comover absorption scenario [ig(l . The different symbols stand for the NA50 data |2(|| from 
the year 2000 (analysis A,B,C) while the dashed histogram is the UrQMD result [ig- 

works see [la, [l3, [3 Backward channels 'charm + anticharm meson — ^ charmonia + meson' are treated via detailed 
balance in a more schematic interaction model with a single parameter or matrix element |Mp that is fixed by the 
J/^ suppression data from the NA50 collaboration at SPS energies (cf. Ref. [iflj). 

We recall that charmonium suppression had been proposed as "the clearest" QGP-signature and the community 
has been riding on this folklore for almost two decades. Hence, these detailed comparisons come as a shock: The 
Ph + Ph results at 160 A-GeV, both from UrQMD and HSD transport calculations are well in line with the data of 
the NA50 Collaboration in Fig. [3 where the J/^ suppression is shown as a function of the transverse energy Et- 
The solid line stands for the HSD result within the comover absorption scenario while the various data points reflect 
the NA50 data from the year 2000 (analysis A,B,C). The data have moved so that they agree now with the HSD and 
UrQMD calculations [ig (dashed histogram in Fig. \^. We mention that there might be alternative explanations 
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FIG. 3: The calculated J/'^ multiplicity per binary collision - multiplied by the branching to dileptons - as a function of 
the number of participating nucleons, N^grt, in comparison to the preliminary data from the PHENIX collaboration j23| for 
Au + Au and pp reactions (taken from Ref. |19|'l. 



for J/^ suppression, as discussed in Refs. |2ll 12^ and/or further (dissociation) mechanisms not considered here. 
However, for the purposes of the present review it is sufficient to point out that 

a) no sign of unusual physics can be related to the J/^ data and 

b) the models employed here (cf. Figs. 6 and 7 in [l3|) use upper limits for the dissociation cross sections and do 
not lead to a sizeable re-creation of charmonia hy D + D channels at SPS energies. 

At RHIC central Au + Au collisions at ^/s = 200 GeV will, however, produce multiplicities of open charm pairs 
about 2 orders of magnitude larger than at 160 A- GeV, such that a much higher J/^ reformation rate (~ A'^g) is 
expected (cf. Ref. [22|). At RHIC top energy, ^/s — 200 GeV, the J/^ comover dissociation may no longer be 
important, since the charmonia dissociated in this channel are approximately recreated in the backward channels. 
Hence, the J/'^ dissociation at RHIC should be less pronounced than at SPS energies. 

The preliminary data of the PHENIX Collaboration 23] allow for a first glance at the situation encountered in 
Au + Au collisions at ^/s — 200 GeV. Fig. 13 shows the J/^ multiplicity per binary collision as a function of the 
number of participating nucleons Npart in comparison to the data at midrapidity. The statistics is quite limited; 
thus no final conclusion can presently be drawn. However, the data neither suggest a dramatic enhancement of J/\I' 
production nor a complete 'melting' of the charmonia in the QGP phase. 

C. Historical Interlude 

Hydrodynamic flow and shock formation has been proposed early |2J, I2ItI | as the key mechanism for the creation of 
hot and dense matter during relativistic heavy-ion collisions. The full three-dimensional hydrodynamical flow problem 
is much more complicated than the one-dimensional Landau model |26[ : the 3-dimensional compression and expansion 
dynamics yields complex triple differential cross-sections, which provide quite accurate spectroscopic handles on the 
equation of state. The bounce-off, vi{pt), the squeeze-out, V2{pt), and the antiflow |23, Isl, |23, ■3(1 , ^1] (third flow 
component j32i. ii33j ) serve as differential barometers for the properties of compressed, dense matter from SIS to RHIC. 
The most employed flow observables are: 

- ^ (14) ■ 

Here, Px denotes the momentum in a;-direction, i.e. the transversal momentum within the reaction plane and Py the 



transversal momentum out of the reaction plane. The total transverse momentum is given as pt — \/Px +Pv' ^^^ 

z-axis is in the beam direction. Thus, flow vi measures the "bounce-off", i.e. the strength of the directed flow in the 
reaction plane, and V2 gives the strength of the second moment of the azimuthal particle emission distribution, the 
so-called " squeeze-out" . jHElElilllimilll. In particular, it has been shown jl^, IHIllilH EH that the 
disappearence or "collapse" of flow is a direct result of a first order phase transition. 

Several hydrodynamic models have been used in the past, starting with the one- fluid ideal hydrodynamic approach. 
It is well known that the latter model predicts far too large flow effects. To obtain a better description of the 
dynamics, viscous fluid models have been developed |2J, 1^3) 1^- In parallel, so-called three-fluid models, which 
distinguish between projectile, target and the fireball fluids, have been considered [33. Here viscosity effects appear 
only between the different fluids, but not inside the individual fluids. The aim is to have at our disposal a reliable, 
three-dimensional, relativistic three-fluid model including viscosity 35 , 36] . 

Flow can be described very elegantly in hydrodynamics. However, also consider microscopic multicomponent (pre-) 
hadron transport theory, e.g. models like qMD 'ss'], IQMD [H, UrQMD [13] or HSD j4lj, as control models for 
viscous hydro and as background models to subtract interesting non-hadronic effects from data. If Hydro with and 
without quark matter EoS, hadronic transport models without quark matter - but with strings - are compared to 
data, can we learn whether quark matter has been formed? What degree of equilibration has been reached? What 
does the equation of state look like? How are the particle properties, self energies, cross sections changed? 

To estimate systematic model uncertainties, the results of the different microscopic transport models also have to be 
carefully compared. The two robust hadron/string based models, HSD and UrQMD, are considered in the following. 
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FIG. 4: Left: Sideward flow p^ = wi • pr of K, A and p's at 6 A-GeV as measured by E895 in semi-central collisions at the 
AGS. Right: The same directed flow data for p and A compared to UrQMD calculations for b < 7 f m |45j . 
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D. Review of AGS and SPS results 



Microscopic (pre-)hadronic transport models describe the formation and distributions of many hadronic particles 
at AGS and SPS rather well |43|. Furthermore, the nuclear equation of state has been extracted by comparing to flow 
data which are described reasonably well up to AGS energies ^32, 43, 44, 45, 46, 47]. Ideal hydro calculations, on the 
other hand, predict far too much flow at these energies |3J|. Thus, viscosity effects have to be taken into account in 
hydrodynamics. 

In particular, ideal hydro calculations are factors of two higher than the measured sideward flow at SIS 34] and AGS, 
while the directed flow Px/iTi measurement of the E895 collaboration shows that the p and A data are reproduced 
reasonably well |45| (Fig. 0J in UrQMD, i.e. in a hadronic transport theory with reasonable cross-sections, i.e. realistic 
mean-free-path of the constituents. 
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FIG. 6: Time evolution of directed flow px/N as a function of rapidity for Au+Au collisions at 8 A-GeV in the one-fluid model. 
Left: Hadronic EoS without phase transition. Right: An EoS including a first order phase transition to the QGP p^ . 



Only ideal hydro calculations predict, however, the appearance of a so-called "third flow component" [S^] or 
"antiflow" jig in central collisions (cf. Fig. (S}. We stress that this only holds if the matter undergoes a first order 
phase transition to the QGP. The signal is that around midrapidity the directed flow, Px{y), of protons develops a 
negative slope! In contrast, a hadronic EoS without QGP phase transition does not yield such an exotic "antiflow" 
(negative slope) wiggle in the proton flow wi(y). The ideal hydrodynamic time evolution of the directed flow, Px/N, 
for the purely hadronic EoS (Fig. Ell.h.s.) does show a clean linear increase oi px{y), just as the microscopic transport 
theory (Fig. ^r.h.s.) and as the data (Fig. ^l.h.s.). For an EoS including a first order phase transition to the 
QGP (Fig. Elr.h.s.) it can be seen, however, that the proton flow vi ^ Px/pr collapses; the collapse occurs around 
midrapidity. 

This observation is explained by an antiflow component of protons, which develops when the expansion from the 
plasma sets in p9l |. 

The ideal hydrodynamic directed proton flow Px (Fig. E)) shows even negative values between 8 and 20 A-GeV. An 
increase back to positive flow is predicted with increasing energy, when the compressed QGP phase is probed. But, 
where is the predicted minimum of the proton flow in the data? The hydro calculations suggest this "softest point 
collapse" is at E^ab ~ 8 A-GeV. This has not been verified by the AGS data! However, a linear extrapolation of the 
AGS data indicates a collapse of the directed proton flow at E^ab ~ 30 A-GcV (Fig. [TJ. 

Recently, substantial support for this prediction has been obtained by the low energy 40 A-GeV SPS data of the 
NA49 collaboration [SJ (cf. Fig. |SJ). These data clearly show the first proton "antifiow" around mid-rapidity, in 
contrast to the AGS data as well as to the UrQMD calculations involving no phase transition (Fig. O. Thus, at 
bombarding energies of 30-40 A-GeV, a first order phase transition to the baryon rich QGP most likely is observed; 
the first order phase transition line is crossed (cf. Fig. 0. This is the energy region where the new FAIR- facility at 
GSI will operate. There are good prospects that the baryon flow collapses and other first order QGP phase transition 
signals can be studied soon at the lowest SPS energies as well as at the RHIC fragmentation region j/ > 4 — 5. These 
experiments will enable a detailed study of the first order phase transition at high /xs and of the properties of the 
baryon rich QGP. 
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FIG. 7: Measured SIS and AGS proton dp^c / dy-slope data compared to a one-fluid hydro calculation. A linear extrapolation of 
the AGS data indicates a collapse of flow at E^at ~ 30 A-GeV, i.e. for the lowest SPS- and the upper FAIR- energies at GSI 
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FIG. 8: VI at SPS, 40 A-GeV and 158 A-GeV |s3] . The proton antiflow is observed in the NA49-experiment even at near 
central collisions, which is in contrast to the UrQMD-model involving no phase transition (Fig. |^. 



II. PROTON ELLIPTIC FLOW COLLAPSE AT 40 A GEV - EVIDENCE FOR A FIRST ORDER PHASE 

TRANSITION AT HIGHEST NET BARYON DENSITIES 



At SIS energies microscopic transport models reproduce the data on the excitation function of the proton elliptic 
flow V2 quite well: A soft, momentum-dependent equation of state [S^, ISa . |5J| seems to account for the data. The 
observed proton flow V2 below ^ 5 A-GeV is smaller than zero, which corresponds to the squeeze-out predicted by 
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FIG. 9: Proton and pion flow vi = Px/pr at 40 A-GeV as obtained within the UrQMD model. No proton antiflow is generated 
in this hadronic transport theory without phase transition. 
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FIG. 10: V2 excitation function of protons at the AGS. The E895-E877 data show the transition from squeeze-out to in-plane 
proton elliptic flow at 4-5 A- GeV; the UrQMD calculations show a strong sensitivity to the EoS |45| . 



hydrodynamics long ago "24", "251, ll^, IH, p?, "s^, "sT] . The AGS data (Fig. ^J) exhibit a transition from squeeze-out to 
in-plane flow in the midrapidity region. The change in sign of the proton f 9 at 4-5 A-GeV is in accord with transport 
calculations (UrQMD calculations in Fig. [TUl 45]; for HSD results see [M^l)- At higher energies, 10-160 A-GeV, a 
smooth increase of the flow V2 is predicted from the hadronic transport simulations. In fact, the 158 A-GeV data of 
the NA49-collaboration suggest that this smooth increase proceeds between AGS and SPS as predicted. Accordingly, 
UrQMD gives considerable ( 3%) V2 flow for midcentral and peripheral protons at 40 A-GeV - Fig. ITTI 45] . 

This is in strong contrast to recent NA49 data at 40 A-GeV (cf. Fig. I12|) : A sudden collapse of the proton flow 
is observed for midcentral as well as for peripheral protons. This collapse of V2 for protons around midrapidity at 
40 A-GeV is very pronounced while it is not observed at 158 A-GeV. The UrQMD calculations, without a phase 
transition, show a robust, but wrong 3% flow of protons - in str ong contrast to the data. 

A dramatic collapse of the flow vi is also observed by NA49 j5l|. again around 40 A-GeV, where the collapse of 
V2 has been observed. This is the highest energy - according to 1, 2] and Fig. n- at which a flrst order phase 
transition can be reached at the central rapidities of relativistic heavy-ion collisions. We, therefore, conclude that 
a first order phase transition at the highest baryon densities accessible in nature has been seen at these energies in 
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FIG. 11: Elliptic flow V2 of protons (lower frame) and pions (upper frame) versus rapidity for Pb-fPb collisions at 40 A-GeV 
from the UrQMD calculations ^ . 
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FIG. 12: Elliptic flow V2 of protons versus rapidity at 40 A-GeV Pb-|-Pb collisions |5H | as measured by NA49 for three centrality 
bins: central (dots), mid-central (squares) and peripheral (triangles). The solid lines are polynomial fits to the data |51ll . 

Pb+Pb collisions. Moreover, Fig. E| shows that the elhptic flow clearly distinguishes between a first order phase 
transition and a crossover. 



III. STRONG COLLECTIVE FLOW AT RHIC SIGNALS A NEW PHASE OF MATTER 



The rapid three-body thermalization found by Xu and Greiner (cf. Sec. IV. A) justifies a posteriori the use 
of hydrodynamical calculations for the time evolution of the complex four-dimensional expansion of the plasma. 
However, there is no justification for the use of simple ideal hydrodynamics (i.e. neglecting the important transport 
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FIG. 13: Time evolution of the momentum anisotropy V2 for a cross over EoS. Both, in- and off-equilibrium calculations (upper 
lines) show no drop of the flow V2. The calculation with a first order phase transition off-equilibrium (lower line) shows - on 
the other hand - the collapse of U2-fiow of protons. Hence, the collapse is only possible for a first order transition |5q| . 



and 200 GeV suggest energy independent V2{ti) distributions. Furthermore, the observed distribution has a triangular 
shape. This finding is in strong disagreement with Bjorken boost invariant hydro predictions |12.l57l|. which fit only 
the midrapidity region. The predicted average proton U2-values obtained from the SPHERIO hydro code with NEXUS 
initial conditons (Fig. I14L jSgl) are by factors of two higher than simple smooth initial state hydrodynamic calculations. 
This indicates that ideal hydro with naive smooth initial conditions - as used by many authors - do not describe but 
rather fit the data. Strong viscosity effects can play a role for particles with pt< 1.2 GeV/c: a decent description 
of the dynamics requires, however, relativistic viscous hydro simulations 35, 361 l59l|. The NexSpherio simulations 
(Fig. ^5 58] ) predict very large event-by-event fiuctuations of V2 caused by the strongly fluctuating initial conditions 
(given by NEXUS). Is this in accord with data? What about the effect on the event plane determination? 
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FIG. 14: Elliptic flow V2 of protons as a function of pr for different centralities. NexSpherio ideal hydro results exhibit about 
a factor of two higher proton-W2-values than the PHENIX data up to pT = 1 GeV/c. Furthermore, the fluctuations in the 
NEXUS initial conditions for the SPHERIO- ideal hydro simulation reflect > 50 % event-by-event fluctuations of the proton-W2 
values. Only large viscosity effects can damp out the too large V2 flow itself as well as the fluctuations in V2- 

Microscopic transport simulations of particle yields, dN/dy distributions, etc. give a good description of the RHIC 
Au-|-Au data [19]. The HSD and UrQMD transport approaches are based on string, quark, diquark {q,q,qq,qq) 
as well as hadronic degrees of freedom. At RHIC, UrQMD and HSD yield reasonable abundances of light hadrons 
composed of u, d, s quarks ^. Do they also predict the collective flow properly? The UrQMD prediction is clearly not 
compatible with the measured 6% elliptic flow - it is sizeably underestimated JG^]. When shortening the formation 
time l60| one can get the model results closer to the data (Fig. I15|l but more additional initial pressure - needed to 
create the missing extra flow - is not justified in the hadronic transport models. At high transverse momenta {pt ~ 2 



^ For a more recent survey on hadron rapidity distributions from 2 to 160 A-GeV in central nucleus-nucleus collisions within the HSD 
and UrQMD transDort aDDroaches we refer the reader to Ref. 14211. 
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FIG. 15: High pr- V2 values as a function of Npart as measured by STAR are compared to HSD calculations. The V2 data are 
more than 5 times higher than the HSD model Dredictions for the most central collisions 1611 . 
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FIG. 16: The elliptic flow V2 for charged hadrons (HSD, solid lines) as a function of pseudorapidity i] (upper part) and as a 
function of the number of 'participating nucleons' Apart for \ri\ < 1 (lower part) for Au + Au collisions at ^/s = 200 GeV, in 
comparison to the preliminary 'hit-based analysis' data of the PHOBOS Collaboration |6^ . Note, that in spite of the shortened 
formation, r ~ 0.8 fm/c, HSD clearly underpredicts the data even at the moderate pr-values dominating the pr integrated 
U2-values shown here. At higher pT-values, Fig. 1151 the discrepancy to the data is more dramatic. 



GeV/c) the W2-fl-Ow is underestimated even by a factor of three (Fig. I15|l in the HSD model '6l!|. We mention that 
the microscopic quark-gluon-string model inserts in addition short distance vector repulsion in order to achieve high 
flow values |6^ . 

Is a transport approach - based on strings and hadronic degrees of freedom - adequate in the initial stage of 
nucleus-nucleus collisions at RHIC energies, where the quark-gluon plasma is formed? Well, the particle abundancies 
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show a rather smooth evolution from SIS to RHIC energies |64ll65| . However, the effective partonic degrees of freedom 
in the initial phase are needed to supply the large pressure obviously needed to describe the elliptic flow at RHIC 
energies. Even 'early' hadron formation - as in HSD with r — 0.8 fm/c ~ and 'large' (pre-)hadronic interaction cross 
sections are insufficient to explain the V2 flow data. This is demonstrated in Fig. E| which shows the calculated 
elliptic flow V2 from HSD for charged hadrons (solid lines) as a function of the pseudorapidity rj (upper part) and as 
a function of the number of 'participating nucleons' Nnnrt (lower part) for \r]\ < 1 in comparison to the preliminary 



'hit-based analysis' data of the PHOBOS Collaboration f63]. The HSD results are very similar to those of the hadronic 
rescattering model by Humanic et al. £^ .BTj and agree with the calculations by Sahu et al. |68l | performed within 
the hadron-string cascade model JAM [63 . 



IV. EARLY THERMALIZATION AT RHIC - EVIDENCE FOR A NEW PHASE 
A. Elastic and inelastic multi-particle collisions in a parton cascade 

To describe the early dynamics of ultrarelativistic heavy-ion collisions and to address the crucial question of ther- 
malization and the early pressure build up at RHIC, unexpectedly high elastic parton cross sections have been assumed 
in parton cascades [T^, Izj in order to reproduce the elliptic flow V2{pt) seen experimentally at RHIC. These cross 
sections are about 1/9 of the baryon-baryon total cross section (^ 45 mb) or 1/6 of the meson-baryon cross section 
(~ 30 mb), such that the effective cross section for the constituent quarks and antiquarks is roughly the same in the 
partonic and hadronic phase, however, tenfold higher than the cross sections calculated in pQCD. 

It has been a great puzzle until recently, when Xu and Greiner developed a consistent three-body kinetic parton 
cascade algorithm JT^jIZSl- These stochastic inelastic ('Bremsstrahlung') 2^3 and 3 — > 2 collisions (gg ^-> ggg) 
drive early thermalization, rather than the two-body elastic collisions, which are too strongly forward peaked. A 
quantitative understanding of the early dynamical stages of deconfined matter is finally in sight. Parton cascade 
analyses incorporating only binary 2 <-> 2 pQCD scattering processes can not build up thermalization and early 
quasi-hydrodynamic behaviour necessary for achieving sufflcient elliptic flow. The importance of inelastic reactions 
was raised in the so called 'bottom-up thermalization' picture |7J|; gluon multiplication leads to a much faster 
equilibration. 



x=0.2 fm/c 
1:= 1.0 fm/c 
i:=3.0 fm/c 
T=5.0 fm/c 




1 2 3 
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FIG. 17: Time evolution of the midrapidity transverse momentum spectrum for a central RHIC collision from the three- 
dimensional Xu-Greiner-parton cascade simulation [73 . rZ^ (prelimenary results are from i75] - including transverse expansion) . 
Only partons residing in a central cylinder of radius 7? < 5 fm are plotted. The initial off-equilibrium conditions are given 
by a minijet distribution with corresponding overlap function in space-time. At t = 0, only minijets with pr > 2 GeV/c are 
populated. Energy degradation to lower momenta proceeds by rapid gluon emission within the first fm/c. (Quasi-)kinetic and 
chemical equilibrium is found up to 4 fm/c. Here longitudinal and transversal hydrodynamic work is at action, resulting in a 
rapid lowering of the temperature by a factor of two. On the other hand, a small fraction of the initial non-equilibrium high 
momentum power-law tail of the mini-jet production survives even in this central cell. 
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Xu and Greiner consider - besides gluon- and quark- two-body elementary parton-parton scatterings - three-body 
processes gg <-> ggg in leading-order pQCD. They employ effective Landau-Pomeranchuk-Migdal suppression and 
standard screening masses for the infrared sector of the scattering amplitude. 

The early stage of gluon production in the Xu-Greiner (X-G) approach (Fig. I17|l leads to a rapid kinetic equilibration 
of the momentum distribution as well as to a rather abrupt lowering of the temperature by soft gluon emission. Detailed 
balance among gain and loss contributions is reached rapidly, too. The later, slower time evolution is then governed 
by chemical equilibration of the quark degrees of freedom. The X-G cascade does allow to study in detail RHIC- 
collisions with various initial conditions like minijets or color glass condensate (CGC). Fig. El depicts a preliminary 
calculation [75j using minijet-initial conditions. 

Thermalization and chemical equilibration - as proposed in the bottom-up scenario |74| - can thoroughly be tested 
within this approach. The impact parameter dependence on the transverse energy is used to understand elliptic and 
transverse flow at RHIC within this new kinetic parton cascade with inelastic 3 — » 2 and 2 — j 3 interactions. 

V. HOW MUCH QUENCHING OF HIGH pr HADRONS IS DUE TO (PRE-) HADRONIC FINAL 

STATE INTERACTIONS? 

A (mini-)jet at RHIC can produce hard particles, with pT above 5 GeV/c, but must also form soft particles with 
Pt around 2 GeV/c. Jets produced in the center of the plasma zone have to pass first through the parton phase at 
very high temperatures, then through the correlated diquark and constituent quarks and finally through the hadronic 
phase that has build up preferentially close to the surface of the fireball. Very high p^ jets with 7 > 10 materialize 
only far outside of the plasma. Most of the jets - observed at RHIC - are at pt ~ 4 — 5 GeV/c. More than 50% of 
the leading jet particles at pt ^ 5 GcV/c are baryons. Pion jets of 5 GeV have a 7 = 35, i.e., they form far outside 
the plasma. However, HSD-PYTHIA-calculations 76] show that many pious stem from decay ing rho-jets. But, p's 
and protons of 5 GeV have 7 = 5. Thus, p and p-jets hadronize with roughly 50% probability |£l|jlZ3 while passing 
through the expanding bulk matter. All partonic and hadronic models have failed by factors of 5-10 to predict the 
observed high baryon abundance. 

The PHENIX JT^I and STAR 79] collaborations reported a suppression of meson spectra for transverse momenta pr 
above ~ 3 GeV/c. This suppression is not observed in d-|-Au interactions at the same bombarding energy per nucleon 
|80l Isil ] and presents clear evidence for the presence of a new form of matter. However, it is not clear at present how 
much of the observed suppression can be attributed to (pre-)hadronic interactions (FSI) [73 . (In-)elastic collisions of 
(pre-)hadronic high momentum states with some of the bulk (pre-)hadrons in the fireball can contribute in particular 
to the attenuation of pr ~ 5 GeV/c transverse momentum hadrons at RHIC 76] : Most of the medium momentum 
(pre-)hadrons from a ±5 GeV/c double jet will materialize inside the dense plasma; their transverse momenta being 
0-4 GeV/c. The particles are dominantly p's, K's and baryons at pT > 2.5 GeV/c - hence their formation time is 
7Ti? « 4 fm/c in the plasma rest frame. The time for color neutralization can also be very small 82] for the leading 
particle due to early gluon emission. 

The (pre-) hadronic interactions with the bulk of the (pre-)hadronic comovers then must have clearly an effect: they, 
too, suppress the p^-spectrum [73. (In)elastic reactions of the fragmented (pre-)hadrons with (pre-)hadrons of the 
bulk system cannot be described by pQCD: The relevant energy scale ^/s is a few GeV. Such (in-)elastic collisions 
are very efficient for energy degradation since many hadrons with lower energies are produced. On the average, 1 to 2 
such interactions can account for up to 50% of the attenuation of high px hadrons at RHIC [73. Hence, the hadronic 
fraction of the jet-attenuation had to be addressed. 

In Ref. |6j the HSD transport approach 4r] is employed. Moderate to high transverse momenta (> 1.5 GeV/c) 
are incorporated by a superposition oi p + p collisions described via PYTHIA [ij]. In Au+Au collisions, the formation 
of secondary hadrons is not only controlled by the formation time Tf , but also by the energy density in the local rest 
frame. In 61], hadrons are not allowed to be formed if the energy density is above 1 GeV/fm''^. The interaction of 
the leading and energetic (pre-)hadrons with the soft hadronic and bulk matter is thus explicitly modeled. 

Figs. 18, 19 show the nuclear modification factor [Tg 

for the most central (5% centrality) Au+Au collisions at RHIC. The Cronin enhancement is visible at all momenta. 
Hadron formation time effects do play a substantial role in the few GeV region, since heavier hadrons (K*'s, p's, 
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FIG. 18: The suppression factor Raa (3) of charged hadrons at 5% (10%) central Au + Au coUisions (y^=200 GeV) at 
midrapidity (hatched band). The experimental data are from Refs. |23 . ISq and show clearly that an additional partonic 
suppression is needed (taken from Ref. |6l|). 



protons) are formed 7 times earlier than the rather light pions in the cms frame at fixed transverse momentum due 
to the lower Lorentz boost 7 < 5. It was shown in i61i| that for transverse momenta pr > 6 GeV/c the interactions of 
formed hadrons are not able to explain the attenuation observed experimentally. However, the ratio -Raa is influenced 
by interactions of formed (pre-) hadrons in the pT = I . ■ -5 GeV /c range 61]. A similar behaviour has also been found 
in UrQMD simulations l83|. 

As pointed out before, the suppression seen in the calculation for larger transverse momentum hadrons is due to 
the interactions of the leading (pre-)hadrons with target/projectile nucleons and the bulk of low momentunr hadrons. 
It is clear that the experimentally observed suppression can not be quantitatively described by the (pre-)hadronic 
attenuation of the leading particles |61,]. The ratio i?AA (3) decreases to a value of about 0.5 at 5 GeV for central 
collisions, whereas the data are around i?AA ~ 0.25. 

To check, how robust this HSD estimate is, alternative models for the leading pre-hadron cross section have been 
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FIG. 19: Same as Fig. 1181 but with a leading cross section according to eq. @ for the perturbative high pt particles. The 
observed attenuation is more than double the value! 
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studied in "61] by adopting a time-dependent, color-transparency- motivated cross section for leading pre-hadrons |84| 

CTload(\/s, t) = Crhad(\/s) (4) 

for T — To < Tf, where tq denotes the actual production time, Tf the formation time, after which the ftiU hadronic 
cross section is reached. 

Within this scenario the attenuation is 35% at pT ^ 5 GeV/c (see Fig. I19f) . while the data show more than 
double the attenuation. Thus, (pre-)hadronic jet interactions cannot provide a quantitative explanation for the jet 
suppression observed. They do provide, however, a sizable (30 — 50%) contribution to the jet quenching. 

For particles observable with momenta pt > 4 GeV/c, the HSD transport calculation predicts that still 1/3 of 
the final observed hadrons have suffered one or more interactions, whereas the other 2/3 escape freely, i.e., without 
any interaction (even for central collisions). This implies that the final high px hadrons originate basically from the 
surface. 

A. Angular Correlations of Jets — Can jets fake the large W2-values observed? 

Fig. 1201 [73 shows the angular correlation of high pT particles {pxNTrig = 4 ... 6 GeV/c, pT = 2 GeV . . .prNTrig, 
\y\ < 0.7) for the 5% most central Au-|-Au collisions at ^/s = 200 GeV (solid line) as well as pp reactions (dashed line) 
from the HSD-model |7Q] in comparison to the data from STAR for pp collisions [83. Gating on high px hadrons 
(in the vacuum) yields 'near-side' correlations in Au-|-Au collisions close to the 'near-side' correlations observed for 
jet fragmentation in the vacuum (pp). This is in a gre ement with the experimental observation 86] . However, for 
the away-side jet correlations, the authors of Ref. [2a| get only a ^50% reduction, similar to HIJING, which has 
only parton quenching and neglects hadron rescattering. Clearly, the observed jSQj complete disappearance of the 
away-side jet (Fig. 1201 cannot be explained in the HSD- (pre-) hadronic cascade even with a small formation time of 
0.8 fm/c. Hence, the correlation data provide another clear proof for the existence of the bulk plasma. 

Although (pre-)hadronic final state interactions yield a sizable (< 50%) contribution to the high pT suppression 
effects observed in Au-|-Au collisions at RHIC, ^ 50% of the jet suppression originates from interactions in the plasma 
phase. The elliptic flow, V2, for high transverse momentum particles is underestimated by at least a factor of 3 in 
the HSD transport calculations |6lj | (cf. Fig. I15|) . The experimentally observed proton excess over pions at transverse 
momenta pr > 2.5 GeV/c cannot be explained within the CGG approach |63|; in fact, the proton yield at high 
Pt > 5 GeV/c is a factor 5-10 too small. We point out that this also holds for partonic jet-quenching models. Further 
experimental data on the suppression of high momentum hadrons from d+Au and Au-|-Au collisions, down to ^/s = 
20 GeV, are desperately needed to separate initial state Cronin effects from final state attenuation and to disentangle 
the role of partons in the colored parton plasma from those of interacting pre-hadrons in the hot and dense fireball. 

Can the attenuation of jets oi pr > 5 GeV/c actually fake the observed W2-values at p^ ~ 2 GeV/c? This question 
comes about since due to fragmentation and rescattering a lot of momentum-degraded hadrons will propagate in the 
hemisphere defined by the jets. However, their momentum dispersion perpendicular to the jet direction is so large 
that it could indeed fake a collective flow that is interpreted as coming from the high pressure early plasma phase. 

On first sight. Fig. 1211 shows that this could indeed be the case: the in-plane V2 correlations are aligned with the 
jet axis, the away-side bump, usually attributed to collective V2 flow (dashed line), could well be rather due to the 
stopped, fragmented and rescattered away-side jet! However, this argument is falsified by the out-of-plane correlations 
(circles in Fig. I21|l . The near-side jet is clearly visible in the valley of the collective flow 112 distribution. Note that 
V2 peaks at if = 7r/2 relative to the jet axis! The away-side jet, on the other hand, has completely vanished in the 
out-of-plane distribution (cf. Fig. I22|l ! 

Where are all the jet fragments gone? Why is there no trace left? Even if the away-side jet fragments completely 
and the fragments get stuck in the plasma, leftovers should be detected at momenta below 2 GeV/c. Hadronic models 
as well as parton cascades will have a hard time to get a quantitative agreement with these exciting data! 

We propose future correlation measurements which can yield spectroscopic information on the plasma. 

1. If the plasma is a colorelectric plasma, experiments will - in spite of strong plasma damping - be able to search 
for wake-riding potential effects. The wake of the leading jet particle can trap comoving companions that move 
through the plasma in the wake pocket with the same speed [pT/m) as the leading particle. This can be 
particular stable for charmed jets due to the deadcone effect as proposed by Kharzeev et al [83, which will 
guarantee little energy loss, i.e. constant velocity of the leading D-meson. The leading D-meson will practically 
have very little momentum degradation in the plasma and therefore the wake potential following the D will be 
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FIG. 21: High pr correlations: in-plane vs. out-of-plane correlations of the probe (jet-l-secondary jet fragments) with the bulk 
(«2 of the plasma at pt > 2GeV/c), prove the existence of the initial plasma state (STAR-coUaboration, preliminary). 



One may measure the sound velocity of the expanding plasma by the emission pattern of the plasma particles 
travelling sideways with respect to the jet axis: The dispersive wave generated by the wake of the jet in the 
plasma yields preferential emission to an angle (relative to the jet axis) which is given by the ratio of the leading 
jet particles' velocity, devided by the sound velocity in the hot dense plasma rest frame. The speed of sound for 
a non-interacting gas of relativistic massless plasma particles is Cs ~ -7= ~ 57% c, while for a plasma with strong 
vector interactions, Cg = c. Hence, the emission angle measurement can yield information of the interactions in 
the plasma. 
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FIG. 22: Illustration of jets traveling through the late hadronic s tag e of the reaction. Only jets from the region close to the 
initial surface can propagate and fragment in the vacuum 123. l76l l89l l9(jl . The other jets will interact with the bulk, resulting 
in wakes with bow waves travelling transversely to the jet axis. 

VI. SUMMARY 

The NA49 collaboration has observed the collapse of both, vi- and W2-collective flow of protons, in Pb+Pb collisions 
at 40 A-GeV, which presents first evidence for a first order phase transition in baryon-rich dense matter. It will be 
possible to study the nature of this transition and the properties of the expected chirally restored and deconfined 
phase both at the forward fragmentation region at RHIC, with upgraded and/or second generation detectors, and 
at the new GSI facility FAIR. According to Lattice QCD results [l|, Q , the first order phase transition occurs for 
chemical potentials above 400 GeV. Fig. ^| shows that the elliptic flow clearly distinguishes between a first order 
phase transition and a crossover. Thus, the observed collapse of flow, as predicted in [24.i25|. is a clear signal for a 
first order phase transition at the highest baryon densities. 

A critical discussion of the use of collective flow as a barometer for the equation of state (EoS) of hot dense matter 
at RHIC showed that hadronic rescattering models can explain < 30% of the observed flow, V2, ioi pr > 2 GeV/c. We 
interpret this as evidence for the production of superdense matter at RHIC with initial pressure way above hadronic 
pressure, p > 1 GeV/fm^. 

The fluctuations in the flow, vi and W2, should be measured. Ideal Hydrodynamics predicts that they are larger 
than 50 % due to initial state fluctuations. The QGP coefficient of viscosity may be determined experimentally from 
the fiuctuations observed. 

The connection of V2 to jet suppression is examined. It is proven experimentally that the collective flow is not faked 
by minijet fragmentation and that the away-side jet suppression can only partially (< 50%) be due to pre-hadronic 
or hadronic rescattering. 

We propose upgrades and second generation experiments at RHIC, which inspect the flrst order phase transition 
in the fragmentation region, i.e. at /is « 400 MeV (y « 4 — 5), where the collapse of the proton flow analogous to 
the 40 A- GeV data should be seen. 

The study of Jet- Wake-riding potentials and Bow shocks caused by jets in the QGP formed at RHIC can give 
further clues on the equation of state and transport coefficients of the Quark Gluon Plasma. 
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